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We study the recent Floquet-realisation of the Harper-Hofstadter model in a gas of cold bosonic 
atoms. We study in detail the scattering processes in this system in the weakly interacting regime due 
to the interplay of particle interactions and the explicit time dependence of the Floquet states that 
lead to band transitions and heating. We focus on the experimentally used parameters and explicitly 
model the transverse confining direction. Based on transition rates computed within the Floquet- 
Fermi golden rule we obtain band population dynamics which are in agreement with the dynamics 
observed in experiment. Finally, we discuss whether and how photon-assisted collisions that may 
be the source heating and band population dynamics might be suppressed in the experimental 
setup by appropriate design of the transverse confining potential. The suppression of such processes 
will become increasingly important as the experiments progress into simulating strongly interacting 
systems in the presence of artificial gauge fields. 


I. INTRODUCTION 

Time-periodic driving provides a versatile tool with 
which to control the behaviour of quantum systems. 
Such systems, described by Floquet theory, can lead 
to forms of long-time dynamics that differ dramatically 
from those in the absence of any drive. In particular, the 
ability to use time-periodically modulated light fields 
to cause neutral atoms to experience artificial gauge 
fields has sparked significant experimental [IHH] and 
theoretical activity in studies of Floquet sys¬ 

tems in recent years. Whereas the long-time dynam¬ 
ics of a single particle in a periodically driven system 
can often be well understood in terms of an effective 
Hamiltonian . the applicability of perturbative 

high-frequency expansions to driven interacting many- 
particle systems is less established [3T] . Generically, un¬ 
der periodic drive, a closed interacting quantum system 
is expected to heat up to an infinite temperature state 
|321484j unless many-body localisation prevents this [31- 
[37]. In open quantum systems, depending on the spe¬ 
cific form of the system-bath coupling, heating to an in¬ 
finite temperature state may be avoided and Gibbs-like 
steady states can emerge |55HiD] . By properly choosing 
the bath density of states and the coupling to the system 
one may obtain steady states that resemble finite tem¬ 
perature states of non-driven systems [4lti44] . In either 
case the heating processes generically present in Flo¬ 
quet systems might make them unsuited to experimen¬ 
tally study interacting and strongly correlated phases. 
Thus, the need to determine the relevant timescales over 
which the heating processes occur becomes immediately 
apparent. 

In a perturbative treatment of interacting Floquet 
systems one finds that inelastic scattering associated 
with the absorption of quanta Hto becomes possible 
which can lead to heating and band transitions jH |35] . 
For heating due to two-particle collisions independent 
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of the specific lattice and driving protocol considered, 
it is found that the rates scale with the square of the 
corresponding interaction strength and scale linearly in 
the density of particles, and in the case of scattering into 
transverse dimensions can be further modified by chang¬ 
ing the transverse density of states. Thus, appropriate 
design of the Hamiltonian might be used to suppress 
these processes, of paramount importance for experi¬ 
ments as deleterious heating and band transitions can 
limit the achievable lifetimes. Depending on the specific 
setup the microscopic processes leading to such energy 
absorption are slightly different. In all cases they can be 
understood as arising from the explicit time-dependence 
of the Floquet-states and the non-commutativity of the 
interaction with the time-evolution operator of the non¬ 
interacting system |45j . If in addition to the perturba¬ 
tive treatment of the interactions, the time-dependence 
is treated approximately, e.g. in a rotating wave ap¬ 
proximation, or the time-evolution is further expanded 
in small parameters, e.g. the hopping on the lattice, 
the heating rates can be associated with a microscopic 
process, e.g. a particle hopping to a lower energy site 
and converting that potential energy into kinetic energy, 
which would be possible in a time-independent system, 
or a genuine Floquet process in which the particles ab¬ 
sorb energy from the driving field and convert that into 
energy of motion |46j . For the system we consider below 
it is the second process which will be of most interest, it 
being unique to Floquet systems. This process will nat¬ 
urally scale with the amplitude of the driving field, the 
rates typically quadratically within perturbation theory 
in the driving fields for small amplitudes, and as a com¬ 
plicated function for stronger driving. 

In this paper, we study heating processes in a Flo¬ 
quet system consisting of a periodically modulated su¬ 
perlattice potential, as used in a recent experimental 
realisation of the Harper-Hofstadter Hamiltonian (T^j. 
We apply the perturbative formalism we recently de¬ 
veloped to describe Floquet systems used to generate 
artificial gauge fields for cold gas systems |3S|. We will 
treat the non-interacting time-dependent Hamiltonian 
exactly within Floquet theory and treat the particle- 
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interactions perturbatively. Focusing on the effects of 
two-particle scattering processes, we compute the heat¬ 
ing processes that arise from weak two-body interac¬ 
tions and corresponding intra- and inter-band transition 
rates, and compare with experimental observations. We 
further study ways in which these processes could be 
suppressed by appropriate design of the Hamiltonian. 


We begin by introducing the model in Sec. |TT] and 
obtain the single particle Floquet spectrum for the ex¬ 
perimentally used parameters in Sec. |II A[ We proceed 
by introducing the Floquet scattering processes com¬ 
puted within the Floquet-Fermi golden rule in Sec. [Tib 


In Sec. IIC we compute the band population dynamics 
based on the Floquet scattering rates for the experi¬ 
mentally used parameters. Finding reasonable agree¬ 
ment with the experimentally observed dynamics we 
then discuss in Sec. |IID| whether these processes can be 
suppressed by appropriate design of the transverse con¬ 
fining potential. We show that a strong suppression can 
be achieved under sufficiently strong confinement, sug¬ 
gesting a possible way to enable experiments to access 
strongly correlated quantum phases without the delete¬ 
rious heating and repopulation dynamics demonstrated 
to be present in this Floquet system at weak confine¬ 
ment. 


II. MODEL 



FIG. 1. (Color online) Resulting band structure En{kx,ky) 
in units of the recoil-energy = h^k^ / {2M) as a func¬ 
tion of momentum kx and ky in the reduced BZ of the non¬ 
interacting model, Eq. 0. for the parameters 14 = dEr, 
Vxi = O.SEr, Vy = lOEr and k = 0.58huj and huj = 0.72Er. 
The dispersion splits into 4 bands of which the middle two 
bands are touching at the borders of the reduced BZ and 
combined into a single superband. 


mental setup of Ref. [TH is given by 

Vix,y,t) = Vst{x,y) + Vosc{x, y,t) 

= 14 sin^(fcr2;) -I- Vxi sin^(fcra;/2) -|- Vy sin^ (kry) 
+ K [sin(7r/4 -|- kyXjT) cos{(j)o + — kryj2) 

— cos(7r/4 + krxl2) sin(0o — ojt — kry/2)\ 

(3) 


We consider bosons described by a field-operator 4' (x) 
loaded into a two-dimensional (2D) optical lattice with 
both time-dependent and static parts with additional 
transverse confinement. The resulting Hamiltonian is 
given by iJ = iJo(0 + 


iFo(i)= J dx^vl/t(x)A4/(x) 

+ ydxvl/t(x) [Vix,y,t) + V^{z)]^{K), 
dx 4'^' (x)4''f (x)4'(x)4'(x), 


Hint = 


( 1 ) 

( 2 ) 


where the first line of equation 0 describes the kinetic 
energy of the atoms moving in the optical lattice and the 
second line gives the time-dependent in-plane optical 
potential V {x, y, t) and the static transverse confining 
potential V^{z) experienced by the atoms. The last 
part, equation 0, gives the contact interaction between 
atoms of strength g. For the weak contact interactions 
considered here, g = iirh'^as/M describes collisions with 
s-wave scattering length Us- 

We consider the in-plane optical lattice V{x,y,t) to 
be a dynamical superlattice[25j. which for the experi¬ 


where the second line describes the static super lattice 
Vstix, y), consisting of the short lattice in x{y) direction 
of strength Vx{Vy) and the long lattice with strength 14; 
that creates the staggering along the ^-direction. The 
third and fourth line describe the time-dependent part 
Vosc{x, y,t) of strength k oscillating with frequency ui 
which is fixed to be on resonance with the energy stag¬ 
gering along the a:-direction. The time-dependent part 
is chosen such that the second and third line respec¬ 
tively restore tunnelling along each other bond in the 
staggered a:-direction, and such that a homogeneous flux 
of tt/2 is realised in the tight-binding description. The 
phase (j)o is not controlled in the experiment and will be 
fixed in the following to be 4>o = 7r/4; we have checked 
that the results are not changed significantly for a dif¬ 
ferent choice. The transverse confinement potential is 
taken to be either an optical lattice cos^(fc^z) 

or a harmonic trap 144 = which is used in 

the experiment. 


A. Single-particle states and dispersion 

We proceed to obtain the single-particle Floquet 
spectrum of the non-interacting time-dependent Hamil¬ 
tonian iFo(0- Due to the separability of the non¬ 
interacting Hamiltonian we may focus on the parti- 
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cle spectrum of the in-plane {x-y) motion, i.e. wave- 
functions take the form ^{x,y,z,t) = {x,y,t)'^^{z) 

and energies are given by En,k,m^ = . As 

the Hamiltonian is invariant under discrete temporal 
and spatial translations the solutions take the form of 
Floquet-Bloch waves. 

We note that the static potential has the following 
symmetries, Vst{x + 2a,y) = Vst{x,y + a) = Vst{x,y) 
with the lattice spacing a = tt jk^. We may choose a unit 
cell of 2 X 2 sites and each Bloch-band will be split into 4. 
However, the time-dependent part has a lower symme¬ 
try, Voscix + 4a, y,t) = Vosc{x,y + 4a,t) = Vosdx, y,t), 
and naively, this would lead to a real space unit cell of 
4x4 sites and to a splitting into 16 bands. However, 
the unit cell can be reduced by rewriting Vi^sc as 

Vosc{x, y, t) = e“*F(a;, y) + A* {x, y) (4) 

with quasi-periodic F(r -|- Rj) = with G = 

(7r/(2a),7r/(2a)). This allows us to perform a unitary 
gauge-transformation in Floquet-space as done in [5S] 
and obtain a Hamiltonian invariant under translations 
by 2 lattice sites in both x- and y-directions. Conse¬ 
quently, we may keep the unit cell consisting of 2 x 2 
sites with each band split into 4 and exactly expand 
the time-dependent problem in the Bloch-states of the 
time-independent problem as they now share the same 
periodicity. 

We project the full time-dependent Schrodinger equa¬ 
tion on the set of the 4 lowest Bloch-bands of 
which are resonantly coupled by the time-dependent op¬ 
tical potential T4sc- Moreover, in the expansion of the 
Floquet-states, 4'e(t) = 

we keep only a finite number of frequency components 
—M < m < M with sufficiently high M to ensure con¬ 
vergence of both the spectrum and wavefunctions. 

A plot of the resulting band structure is shown in 
Fig. 0 which shows the 4 Harper-Hofstadter bands in 
the reduced Brillouin zone corresponding to the 2x2 
real-space unit-cell. Following the convention in we 
refer to the middle two bands as a single band as they 
are not separated by an energy gap. 


B. Floquet-Scattering 

We treat the effects of the collisions of the atoms given 
by the interaction Hamiltonian 0 in the framework 
of scattering theory for Floquet states, [35j and refer¬ 
ences therein. We employ the Floquet-Fermi golden rule 
(FFGR) to compute transition rates from an initial state 
4'i to final state given by 

= E I I 

Am 

(5) 

where is the time-periodic Floquet mode associ¬ 
ated with double 

brackets denote the extended scalar product given by 


1 

10"^ 


10"^ 
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FIG. 2. (Color online) Overlap matrix element /Am = 
y) | Hmt \ j/)))|^ of the wavefunctions of the 

2 dimensional model for a transition (1,1) —> (3, 3) summed 
over the final state momenta and averaged over the initial 
state momenta as a function of the photon transfer Am nor¬ 
malised to the element at Am = 0. 



(($f I $i)) = 1/T I $i(t)). The sum over Am 

allows the absorption or emission of Am energy quanta 
hw during the scattering process. 

Using the FFGR will preclude the discussion of 
strongly correlated many-body phases, but is appropri¬ 
ate for the regime of weakly-interacting particles con¬ 
sidered here and sufficient to explore the two-particle 
physics that are a potential source of heating and band 
repopulation. 

We apply FFGR Eq. 0 to an initial state contain¬ 
ing particles in bands rii and n 2 with in-plane crystal- 
momentum ki and ^2 of the 2D bandstructure in states 
mf and m| of the transverse potential, i.e. 




at 


4' 


at 


ni,fei,m^ n2,fc2,m| 


I vac) 


( 6 ) 


scattering into a final state containing two particles in 
bands ns and with crystal-momentum ks and /c 4 in 
states m| and m| 


'kf = 'k 


at 

n3,fc3,m 


3 


at 

n4,/c4,m 


e I vac). 

4 ' 


(7) 


In the calculations of the population dynamics in 
Sec. H G we only keep transitions with emission and 
absorption of up to a single energy quantum huj, i.e. 
Am = 0,±1 in the FFGR Eq. 0, as transitions with 
higher energy transfers are successively suppressed as 
shown in Fig. 0. Note that we still have to keep a high 
nnmber M of Floquet modes in the expansion of the 
Floquet states used to compute the rates. The exponen¬ 
tial decay of these transition-matrix elements also jus¬ 
tifies our restriction to the lowest 4 Harper-Hofstadter 
bands. Whereas transitions into the higher bands are 
allowed within FFGR, for the experimental parameters 
they correspond to a transition with Am « II and are 
thus negligible on the time-scales we are interested in. 
This argument applies rather generally and does not rely 
on the precise form of the wavefunctions, but only on 
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the fact that the higher bands are gapped in the static 
Hamiltonian and that the Floquet states are exponen¬ 
tially localised in the frequency domain, thus, leading 
to the observed exponential suppression of the matrix 
elements. 

For the general discussion of the scattering rates in 
Sec. II D| we assume an initial state in the unique low¬ 
est energy state of the transverse direction, i.e. mf = 
= 0 for harmonic confinement, and define a total di¬ 
mensionless rate for transitions ^(ni,n 2 )^{n 3 ,n 4 ) taking 
the initial state to be (ni, fci, 0 ; n 2 , k 2 , 0 ), summing over 
the final state momenta ks and k^ and transverse state 
quantum numbers and m| and averaging the ini¬ 
tial state momenta ki and k 2 over the reduced Brillouin 
zone (BZ) 


r(ni,n2)->(n3,'n4)7 



(iki(ik2(ik3C?k4 


BZ 


( 8 ) 


where we factored out the dimensionful prefactor 
7 o. The prefactor 70 takes the form 70 = 

hET L L hu’’c harmonic confinement and 

T'o “ hE~ ^r LLL for Confinement by an optical 

lattice. 


C. Population dynamics 

We compute the resulting band population dynamics 
based on the Floquet transition rates, Eqn. ([^, for the 
experimental parameters, i.e. Vx = 6Er, Vxi = O.SEr, 
Vy = lOEr, K = O.bSHuj, fvjj = 0 . 72 Er and Er = 
with kr = 7 r/a and a = 0.5 x 767 nm, and for the exper¬ 
imentally used harmonic confinement = ^muioscz'^ 
with Wosc = 12Hz. We use the value of the scatter¬ 
ing length for ®^Rb and the experimental density for 
which N(P/[LxLy) « 20. We assume the particles 
to be homogeneously spread over the BZ and average 
all rates over the initial state momenta as discussed 
above. This allows us to formulate the problem only 
in terms of the band-populations of the 2D Harper- 
Hofstadter model and the transverse direction. The 
resulting dynamics for the fractional band occupations 
riy, = /Ntotai Summed over the states of the 

transverse confining harmonic potential are shown in 
Figi for initial fractional occupations of ni = 0 . 6 , 
712 = 0.3 and 713 = 0.1, all assumed to be initially in 
the = 0 state. The behaviour compares favourably 
with the experimental results, see Fig. 3 in [12]. 

We observe two time scales for the population dy¬ 
namics. On a fast time-scale elastic (Am = 0) collisions 
of the type (1,3) —?> (2,2) and the reverse process re¬ 
distribute particles between the bands while conserving 
the in-plane motion energy, i.e. the transverse motion 
quantum numbers m| = m^ are conserved. As expected 




FIG. 3. (Color online) Dynamics of the fractional band pop¬ 
ulations Uy = /Wotai in the Floqnet realization 

of the Harper-Hofstadter model for the experimental param¬ 
eters with a transverse harmonic confinement summed over 
the mE oscillator quantum number. Colours correspond to 
the bands shown in Fig. 


from Fig.j^these occur considerably faster than the Flo¬ 
quet Am ^ 0 transitions. This initally leads to an in¬ 
crease of the population in band 2 and the corresponding 
decrease of populations in band 1 and 3, while keeping 
all particles in the m^ = 0 states of the transverse poten¬ 
tial. For the initial conditions, these processes achieve 
a quasi-equilibrium on a very short time-scale stopping 
further redistribution. For longer times Floquet transi¬ 
tions (Am = ±1), changing both the Hofstadter bands 
and allowing transfer into higher energy states of the 
tranverse direction mf ^ m^, become relevant. The in¬ 
terplay of the elastic and the inelastic collisions is seen 
to reverse the initial decrease in population of band 3. 
As the higher states of the transverse potential are oc¬ 
cupied the reverse scattering processes become impor¬ 
tant and slow down the population dynamics after a 
time scale of about 10 ms. The further dynamics is 
seen to be slower explained by the fact that the higher 
states of the transverse potential are less strongly inter¬ 
acting. Finally, the system approaches a state in which 
all Harper-Hofstadter bands are equally populated over 
a timescale of 200 ms. 

We remark that even though the Harper-Hofstadter 
band populations approach a steady state, the system 
slowly continues to heat up in the transverse direction. 
The total absorbed energy AE{t) = E{t) — E{t = 0) 
per particle in units of Er is shown in Fig. After an 
initial period of fast energy absorption the rate flattens 
out. This is due to the fact that the interactions be¬ 
come weaker as the particles are transferred to higher 
m^ states. In the beginning all particles occupy m^ = 0 
states and are strongly interacting. The Am = ±1 pho¬ 
ton transitions then establish an equilibrium between 
states at m^ = 0 and m^ = mf, = Huj/Eosc- Conse¬ 
quently, particles in the mf, states would be scattered 
to states in the m^ = 2mf state. The timescale for 
this to happen is controlled by the square of the ratio 
of the matrix elements I Hint \ ) 
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FIG. 4. (Color online) Total absorbed energy AE{t) — 
E{t) — E{t = 0) per particle in units of Er in the Floquet 
realization of the Harper-Hofstadter model for the experi¬ 
mental parameters with a transverse harmonic confinement. 


and I Hint \ /2^mj,/2) which is approximately 

0.06, yielding transition rates 3 orders of magnitude 
lower. We observe that the total absorbed energy per 
particle remains small compared to the bandgaps of the 
undriven system, thus, being consistent with only con¬ 
sidering the 4 lowest bands for these time scales. 

We conclude that the strong inter band population 
dynamics and the associated heating stemming from 
Am = ±1 transitions would make the realization of 
strongly interacting phases exceedingly difficult unless 
proper care is taken to suppress the corresponding pro¬ 
cesses. 


D. Stability and scattering rates 

Based on the conservation of quasi-energy in the 
FFGR Eq. ([^ and the specific experimental parameters, 
one can readily envisage ways in which the transverse 
density of states can be adapted in order to make the 
system stable to the two-particle scattering processes 
with Am = ±1 which cause band transfer and lead to 
heating. 

We consider the cases in which the transverse poten¬ 
tial in the ^-direction is either an optical lattice = 
VzCOS^{k^z) or a harmonic trap = ^moOascZ^■ 

For both cases we investigate whether particles in the 
ground state band[17] of the Harper-Hofstadter model 
can possibly undergo transitions with Am = ±1 de¬ 
pending on the strength of the transverse potential. The 
resulting stability diagram in the case of the optical lat¬ 
tice is shown in Fig. ([^ and in the case of the har¬ 
monic trap in Fig. Q. This stability is based purely 
on kinematic constraints of the scattering process and 
is therefore valid for two-particle scattering beyond the 
applicability of the FFGR. We conclude that by suffi¬ 
ciently strong transverse potential the inelastic 2 par¬ 
ticle single-photon-transfer scattering processes can be 



FIG. 5. (Color online) In (a) stability of the groundstate 
band of the Harper-Hofstadter model to different inelastic 
transitions as indicated in the figure in the case of trans¬ 
verse confinement by an optical lattice = VzCos^{k^z) 
of strength Vz/E^ with E^ = Er- A jump from 0 to 1 in¬ 
dicates the change of allowed to forbidden for the respective 
transition. In (b) the critical confinement strength Vz jE^ 
required to suppress the Am = ±1 transitions as a function 

oiEljEr. 


completely suppressed. 

We begin by discussing the case in which the trans¬ 
verse potential is an optical lattice. For vanishing depth, 
14 = 0, we recover the situation of free particles where 
the system can always absorb energy and is therefore un¬ 
stable. As we increase the lattice depth 14 tbe transition 
with the smallest energy transfer, the intraband scat¬ 
tering transition (1,1) —>■ (1,1), is suppressed first, and 
those with higher energy transfer are consecutively sup¬ 
pressed until finally the last transition, (1,1) —)■ (3,3), 
becomes energetically forbidden and the system is sta¬ 
ble to single photon (Am = ±1) transitions. This hap¬ 
pens when the density of states in the transverse di¬ 
rection vanishes for the required energy transfer, i.e. 
for these parameters when the bandwidth of the low¬ 
est band in the transverse direction becomes smaller 
than huj — 2 [max(E 3 ) — min(iili)]. Additionally, we may 
change the ratio E^fEr which determines the critical 
value of the confinement strength as shown in the right 
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EoscI Er 



EqscI Er 


FIG. 6. (Color online) Scattering rates for 

the band indicated in the figure in the case of transverse 
confinement by an optical lattice = VzCos^{k^z) for 
= Er in panel (a) and E^ = 0.36i?r in (b) corresponding 
to the minimum of observed in Fig. 


FIG. 7. (Color online) Stability of the groundstate band of 
the Harper-Hofstadter model to different inelastic transitions 
as indicated in the figure in the case of transverse confine¬ 
ment by a harmonic potential Vole = ^muJoaeZ^ in panel (a) 
and corresponding scattering rates in panel (b) both as a 
function of the oscillator energy Eoac = fuuJoac in terms of the 
recoil energy Er- 


panel of Fig. ([^ which is minimised at 14 = for 

E^/Er = 0.36. 

The corresponding transition rates for E^ = Er are 
shown in Fig I®. These start from the finite non-zero 
value at 14 = 0 expected for free particles, and increase 
with increasing 14 to reach a maximum. This maxi¬ 
mum occurs at the point at which there is a closing of 
the channel for scattering into higher transverse modes, 
and arises from the characteristic singular density of 
states of the transverse lattice at the edges of the BZ. 
Thus, we see that the transverse optical lattice at first 
increases the scattering rates until it finally completely 
suppresses two-particle single-photon scattering for suf¬ 
ficiently deep lattices. In addition to the lattice depth 
Vz/E^ which determines the form of the density of states 
in the transverse direction, changing the ratio E^/Er 
which effectively rescales the energy axis of the trans¬ 
verse direction compared to the 2D Harper-Hofstadter 
model can be used to control the scattering rates. This 
effect can be understood by the fact that at fixed huj 
reducing E^ requires the transitions to couple to a fi¬ 
nal state containing particles in higher Bloch bands of 
the transverse lattice which have smaller overlap with 


the ground state Bloch band which is assumed to be 
the initial state. This leads to the decrease of the peak 
scattering rate as seen in panels (a) and (b) in Fig. ® 

For the harmonic trap the behaviour is slightly more 
complicated as seen in Fig (®. Again for no transverse 
potential, E^sc = ^osc = 0, we recover the behaviour of 
the free system. As we increase the confinement energy 
Ease transitions become possible and forbidden when¬ 
ever the required transition energy matches an integer 
number of the oscillator energy. For weak confinement, 
this typically includes a range of possible final oscilla¬ 
tor states for a given band-transition which shrinks with 
increasing confinement. The system is stable to single¬ 
photon (Am = ±1) transitions when the energy balance 
equation 2nhuJosc = Eo — Atrans has no solution with in¬ 
teger n for all band transitions. Here we obtain two 
stability regions, one in the intermediate regime around 
ujosc ~ w/S and one for strong confinement Wosc ^ w/2. 

In both cases, by sufficiently strong modifications 
of the transverse density of states, the single photon 
(Am = ±1) transitions can be completely suppressed. 
However, in that regime the next to leading order 
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processes of transitions with Am = ±2 may still be 
present. Due to the localisation of the Floquet modes 
observed in the matrix elements seen in Fig. these 
rates are at least an order of magnitude smaller than 
the Am = ±1 rates. In addition they are suppressed 
by the interaction matrix element of the wavefunc- 
tions of the transverse direction for large energy trans¬ 
fer, e.g. in the case of harmonic confinement we have 

I I I H.nt I « 0.06, 

yielding 3 orders of magnitude compared to the elastic 
(Am = 0) and a factor of 4 compared to the Am = ±1 
rates. In particular, the elastic scattering rates will then 
be at least two orders of magnitude higher than the first 
non-vanishing Floquet rates. 

III. CONCLUSIONS 

We have studied the two-particle scattering processes 
occurring in the Floquet-realisation of the Harper- 
Hofstadter Hamiltonian focusing on the experimental 
setup used in |I2) involving a dynamically modulated 
superlattice potential. Using the Floquet-Fermi golden 
rule we compute scattering rates due to particle interac¬ 
tions. Based on these Floquet transition rates we obtain 
the resulting band population dynamics and compare to 
the experimental results of [12] . The agreement between 
the experimental result and the predictions of this rate 
model suggests that two-particle processes and the Flo¬ 
quet scattering rates might play an important role in 
establishing the band dynamics and the heating rates 
in this system. Of particular relevance to the stabil¬ 
ity of a Bose-Einstein condensate (BEG) in one of the 
Harper-Hofstadter bands, we observe that the scattering 
rates do not depend strongly on the final state momenta. 
Consequently, an initial BEG quickly spreads over the 
BZ and into other subbands. 

We remark that in a different experimental realisation 
of the Harper-Hofstadter model |3S] heating processes 
have been investigated very carefully and two-particle 


scattering was found not to be the limiting factor for 
the lifetimes in that setup. 

More generally, our study provides insight into the 
heating dynamics of a closed quantum system with an 
unbounded dispersion subject to periodic driving. We 
find a timescale over which the system approaches an 
infinite temperature state for the bounded degrees of 
freedom of the in-plane motion and a generically differ¬ 
ent time-scale over which the system then continues to 
heat up in the transverse direction. 

Having established these processes in the Floquet sys¬ 
tem with the geometry of ReflT^ we discussed how these 
rates can be influenced by additional transverse confine¬ 
ment. We studied transverse confinement by both an 
optical lattice and a harmonic potential and concluded 
that by choosing a sufficiently strong confinement in¬ 
elastic single photon transitions with the absorption or 
emission of hw can be suppressed completely. More¬ 
over, in the case of confinement by an optical lattice, 
the scattering rates can be further controlled by choos¬ 
ing an optimal values for the ratio in addition to 

the lattice depth V^. By suppressing the single-photon 
Am = ±1 rates one can achieve a regime in which the 
next order Floquet processes Am = ±2 are at least two 
orders of magnitude smaller than the elastic Am = 0 
rates which are responsible for establishing the strongly 
correlated behaviour of the quantum system. This pos¬ 
sibility to strongly suppress the two-particle inelastic 
scattering rates provides one possible route towards the 
design of future experiments aiming to access strongly 
interacting regimes without deleterious scattering and 
heating. 
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